Astronomy & Astrophysics manuscript no. paper 


© ESQ 2008 


February 2, 2008 





On the systematics of asteroseismological mass determinations of 

PG1159 stars 

L. G. Althaus 1 ' 2 *, A. H. Corsico 1 ' 2 *, S. O. Kepler 3 , and M. M. Miller Bertolami 1 ' 2 ** 



o 
o 

(N 
•*-> 
o 

O 



43 

6 



> 

m 
m 



o 



>< 



1 Facultad de Ciencias Astronomicas y Geofisicas, Universidad Nacional de La Plata, Paseo del Bosque S/N, (1900) La 
Plata, Argentina. 

2 Instituto de Astrofi'sica La Plata, IALP, CONICET-UNLP 

3 Instituto de Ffsica, Universidade Federal do Rio Grande do Sul, 91501-970 Porto Alegre, RS, Brazil 
e-mail: acorsico .althaus ,mmiller@f caglp .unlp. edu. ar ; kepler@if . uf rgs .br 

Received; accepted 

ABSTRACT 

Aims. We analyze systematics in the asteroseismological mass determination methods in pulsating PG 1159 stars. 
Methods. We compare the seismic masses resulting from the comparison of the observed mean period spacings with the 
usually adopted asymptotic period spacings, Allf , and the average of the computed period spacings, AIL. Computations 
are based on full PG1159 evolutionary models with stellar masses ranging from 0.530 to 0.741Mq that take into account 
the complete evolution of progenitor stars. 

Results. We conclude that asteroseismology is a precise and powerful technique that determines the masses to a high 
internal accuracy, but it depends on the adopted mass determination method. In particular, we find that in the case of 
pulsating PG 1159 stars characterized by short pulsation periods, like PG 2131+066 and PG 0122+200, the employment 
of the asymptotic period spacings overestimates the stellar mass by about 0.06 Mq as compared with inferences from 
the average of the period spacings. In this case, the discrepancy between asteroseismological and spectroscopical masses 
is markedly reduced when use is made of the mean period spacing AIL instead of the asymptotic period spacing AII|. 
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1. Introduction 

Pulsating PG1159 stars (or GW Virginis) are evolved hot 
stars that can pose constraints to the stellar evolution the- 
ory of post- Asymptotic Giant Branch (AGB). These vari- 
able stars belong to the population of hydrogen-deficient 
objects characterized by surface layers rich in helium, car- 
bon and oxygen (Werner & Herwig 2006) that are con- 
sidered the evolutionary link between post-AGB stars and 
most of the hydrogen-deficient white dwarfs. The origin of 
most PG 1159 stars is traced back to the occurrence of 
post-AGB thermal pulses: a born-again episode induced 
either by a very late thermal pulse (VLTP) experienced 
by a hot hydrogen-rich white dwarf during its early cooling 
phase — see Herwig et al. (1999), Blocker (2001), Lawlor & 
MacDonald (2003), Althaus et al. (2005), Miller Bertolami 
et al. (2006), or a late thermal pulse (LTP) during which 
hydrogen deficiency is a result of a dredge-up episode, 
(see Blocker 2001). During the VLTP, the convection zone 
driven by the helium flash reaches the hydrogen-rich enve- 
lope of the star, with the result that most of the hydrogen 
content is burnt. 

About a third of spectroscopic PG 1159 stars exhibit 
multipcriodic luminosity variations with periods in the 
range 300 — 3000 s, attributable to global nonradial g- 
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modes pulsation (e.g. Quirion et al. 2007). The presence 
of a pulsational pattern in many PG 1159 stars has allowed 
researchers to infer structural parameters — particularly 
the stellar mass — and the evolutionary status of individ- 
ual pulsators — e.g. Kawaler & Bradley (1994), Kawaler 
et al. (1995), O'Brien et al. (1998), Vauclair et al. (2002) 
and more recently Corsico & Althaus (2006). Stellar masses 
of PG1159 stars can independently be assessed by compar- 
ing the values of log g and log T c g, as inferred from de- 
tailed non-LTE model atmospheres (Werner et al. 1991), 
with tracks coming from stellar evolution modeling, i. e. 
the spectroscopic mass (Dreizler & Heber 1998, Werner & 
Herwig 2006). These two different approaches enable us to 
compare the derived stellar masses. 

Recently, considerable observational and theoretical ef- 
fort has been devoted to the study of some pulsating 
PG 1159 stars. Particularly noteworthy is the work of Fu 
et al. (2007) who have detected a total of 23 frequen- 
cies in PG 0122+200 and Costa et al. (2007) who have 
enlarged to 198 the total number of pulsation modes in 
PG 1159—035, making it the star with the largest number 
of modes detected besides the Sun. Parallel to these obser- 
vational breakthroughs, substantial progress in the theoret- 
ical modeling of PG 1159 stars has been possible (Herwig et 
al. 1999, Althaus et al. 2005, Lawlor & Mac Donald 2006). 
In this sense, the new generation of PG 1159 evolutionary 
models recently developed by Miller Bertolami & Althaus 
(2006) (hereinafter MA06) has proved to be valuable at de- 
riving structural parameters of pulsating PG 1159 on the 
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basis of individual period fits — see Corsico et al. 2007a 
and Corsico et al. 2007b, respectively, for an application to 
the hot pulsating RX J2117. 1+3412 and the coolest mem- 
ber of the class, PG 0122+200. These evolutionary models 
are derived from the complete evolutionary history of pro- 
genitor stars with different stellar masses and an elaborate 
treatment of the mixing and extramixing processes during 
the core helium burning and born again phases. The suc- 
cess of these models at explaining both the spread in surface 
chemical composition observed in PG 1159 stars and the lo- 
cation of the GW Vir instability strip in the log T e g — log g 
plane (Corsico et al. 2006) renders reliability to the infer- 
ences drawn from individual pulsating PG 1159. 

As shown in MA06 the employment of detailed PG 1159 
evolutionary models yields spectroscopical masses that are 
systematically lower — by about 0.05 M — than those 
derived from hydrogen-rich post-AGB tracks (Werner & 
Herwig 2006). Most importantly, the resulting astero- 
seismological masses (as inferred from the period spac- 
ings) are usually 10% higher than the new spectroscop- 
ical masses, except for the hot pulsating PG1159 star 
RX J2117. 1+3412, the spectroscopical mass of which is 
more than 20% higher than the asteroseismological one 
(Corsico et al. 2007a). The mass discrepancy is a clear in- 
dication of the uncertainties weighting upon the mass de- 
termination methods, even though the spectroscopic uncer- 
tainties are of that order. 

In an attempt to understand the persisting discrep- 
ancy between the asteroseismological and spectroscopical 
masses, Miller Bertolami & Althaus (2007) have recently 
shown that previous evolution is not the dominant factor in 
shaping hydrogen-deficient post-VLTP tracks. They con- 
clude that the MA06 PG 1159 tracks are robust enough 
as to be used for spectroscopical mass determinations of 
PG 1159-type stars, unless opacities in the intershell re- 
gion are strongly subestimated. Their results make clear 
that the systematic discrepancy between asteroseismolog- 
ical and spectroscopical masses should not be attributed 
to uncertainties in post-AGB tracks; rather, they call for 
the need of an analysis of possible systematics in the as- 
teroseismological mass determination methods. This is pre- 
cisely the core feature of the present work. Specifically, we 
will concentrate on the usually adopted asymptotic period 
spacing approach (Kawaler et al. 1995; O'Brien et al. 1998; 
Vauclair et al. 2002; Fu et al. 2007) used in most mass 
determinations of individual pulsating PG 1159 stars. The 
advantage of this approach lies in the fact that the mean pe- 
riod spacing of PG1159 pulsators depends primarily on the 
stellar mass (Kawaler & Bradley 1994; Corsico & Althaus 
2006). However, the derivation of the stellar mass using the 
asymptotic predictions may not be entirely reliable because 
they are strictly valid for chemically homogeneous stellar 
models, while PG1159 stars are expected to be chemically 
stratified with strong chemical gradients built up during the 
progenitor star life. We will show that this approach overes- 
timates the seismic mass for those pulsating PG 1159 stars 
on the white dwarf cooling track. We will also show that 
the discrepancy between asteroseismological and spectro- 
scopic masses is markedly alleviated if the average of the 
computed period spacings, instead of the asymptotic ones, 
is used. In the next Section, we summarize the seismological 
tools to infer the stellar mass from the observed mean pe- 
riod spacings. We also describe the evolutionary sequences 
employed. In Sect. 3 and[¥]we present our results and com- 



pare them with other mass determinations methods. We 
close the paper in Sect. [5] by summarizing our findings. 

2. Numerical tools 

The most widely used approach to infer the seismological 
mass of pulsating PG 1159 stars lies in the asymptotic pre- 
dictions of the non-radial pulsation theory (with the no- 
table exception of Kawaler & Bradley 1994 and Corsico et 
al. 2007ab). In the asymptotic limit of very high radial or- 
der k (k 1, i. e., long periods), the g-mode periods of a 
chemically homogeneous stellar model for a given degree I 
and consecutive k are separated by a constant period spac- 
ing AIT| S iven b y (Tassoul et al. 1990) 
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where g is the local gravity and Ti the first adiabatic expo- 
nent (see Hansen & Kawaler 1994). Note that the term in 
brackets is the difference between the real and the adiabatic 
density gradients, that determines buoyancy. The integral is 
taken over the g-mode propagation region. Note that An^ 
is a function of the structural properties of the star via the 
Brunt- Vaisala frequency. The seismological stellar mass is 
constrained by directly comparing An^ as computed from 
Eq. [1] with the observed mean period spacing if the effec- 
tive temperature of the target star is known (by means of 
spectroscopy). Full advantage is taken of the fact that the 
g-mode period spacing of PG 1159 pulsators is mostly sen- 
sitive to the stellar mass and only weakly dependent on the 
stellar luminosity and helium-rich envelope mass fraction 
(Kawaler & Bradley 1994). This feature together with the 
fact that no detailed pulsation calculations are required to 
compute An^ turns the asymptotic period spacing into a 
practical tool to infer the stellar mass of pulsating PG 1159 
stars. 

As mentioned, the asymptotic formula given by Eq. [1] 
is strictly valid for chemically homogeneous stellar models 
and in the limit of high k, i.e. long periods. However, ac- 
cording to the current theory of stellar evolution, PG 1159 
stars are expected to be chemically stratified characterized 
by strong chemical transitions built up during the progeni- 
tor star life. This is illustrated by Fig. [1] which displays the 
inner chemical abundance distribution in a typical PG 1159 
star. Two main chemical transitions, emphasized with gray, 
are easily recognized: an inner C/O interface left by the 
extra mixing episodes that ocurred during central helium 
burning (see Straniero et al. 2003) and an He/C/O interface 
that separates the helium-rich envelope from the carbon- 
oxygen core — modeled by nuclear processing in prior AGB 
and post-AGB stages. Such chemical interfaces produce 
clear and distinctive signatures in N, which are critical for 
the mode-trapping properties of the models. These mode 
trapping features strongly disturb the structure of the pe- 
riod spectrum, thus causing the computed g-mode period 
spacing (Hfc + i — n^) to appreciably depart from uniformity 
(see Kawaler & Bradley 1994 and more recently Corsico & 
Althaus 2005, 2006). 



L. G. Althaus et al.: On the systematics of asteroseismological mass determinations of PG1159 stars 

-i | i | i | i | i | i p 

AIlj (high order k) 

Al I. (low order k) 

An," 



X 




Log (l-MMJ 

Fig. 1. The inner chemical abundance distribution corre- 
sponding to a 0.589 Mq PG 1159 model at log T eB = 5.18. 
The approximate locations of the C/O and He/C/O chem- 
ical transition regions are emphasized with gray. 

A more realistic approach to infer the stellar mass of 
PG 1159 stars that does not suffer from the above men- 
tioned shortcomings is to compare the observed period 
spacing with the average of the computed period spacings, 
AIL; . This quantity is assessed by averaging the computed 
forward period spacings in the same range as the observed 
periods, that is 



AIL 



i^An fc = i^(n fe+1 -n fc ) 



(3) 



where n means the number of observed modes (with m=0) 
of the star. In contrast with the asymptotic approach, the 
assessment of the asteroseismological mass via AIL involves 
the computation of the full adiabatic period spectrum. 
Accurate values of the adiabatic pulsation periods of pulsat- 
ing PG 1159 stars requires the employment of full PG 1159 
evolutionary models that reflect the thermal structure of 
their progenitors (Kawaler et al. 1985). In this work, we em- 
ploy the evolutionary models recently developed by Althaus 
et al. (2005), MA06, and Corsico et al. (2006, 2007b) who 
computed the complete evolution of model star sequences 
with initial masses on the ZAMS (assuming a metallicity 
of Z = 0.02) in the range 1 — 3.75 Mq. These authors 
have followed all of the sequences through the thermally 
pulsing and mass-loss phases on the AGB to the PG 1159 
regime. The evolutionary stages corresponding to the com- 
plete burning of protons shortly after the occurrence of the 
VLTP and the ensuing born-again episode that give rise 
to the H-deficient, He-, C- and O-rich composition char- 
acteristic of PG1159 stars have been carefully followed for 
each sequence. The masses of the resulting remnants span 
the range 0.530 — 0.741 M©. For these PG 1159 evolution- 
ary sequences we have computed £ — 1, g-mode adiabatic 
pulsation periods with the same numerical code and meth- 
ods employed in those works (see Corsico & Althaus 2006 
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Fig. 2. The dipole (£ — 1) asymptotic period spacing (AII|, 
solid line) is compared with the average of period spacing, 
AIT, as a function of the effective temperature for the .53, 
0.589 and 0.741 M Q evolutionary sequences. For AIL; we 
consider short and long periods, i.e., low and high k values 
(dotted and dashed lines, respectively). Stages before and 
after the models reach their highest effective temperature 
are shown. Arrows indicate the direction of evolution. 



for details). In what follows, we will use these evolutionary 
models to compute both the mean AIT and asymptotic 

Am. 



3. Discrepancy between the asymptotic and the 
average of period spacings 

Here we employ the evolutionary models described previ- 
ously to assess the asymptotic period spacing, Allf, and 
the average of the computed period spacings, AIT as given 
by Eqs.[T]and[3J respectively. In Fig. [2J which summarizes 
the main result of our work, we show the run of these two 
quantities for £ — 1 modes (as most detected periodicities 
are triplets) in terms of the effective temperature for se- 
lected stellar masses. To assess the dependence of AIT on 
the period range where the average of the period spacing is 
done, we compute AIT for intervals of short and long peri- 
ods (300-600 s and 900-1500 s, respectively). Different stars 
have different ranges. The resulting AIT in each case are 
denoted by dotted and dashed lines. Note that both AII| 
and AIT decrease as the stellar mass increases. Note also 
that, when the star evolves along the white dwarf cooling 
track, the period spacings increase with decreasing effective 
temperature. This is due to the increasing degeneracy in the 
core as the star cools, causing the Brunt- Vaisala frequency 
to gradually decrease, and the consequent slow increment 
in the periods. 

Most importantly, note from Fig.[2]that, generally, AIT 
turns out to be smaller than AITf. Note also the marked 
dependence of AIT on the period interval where it is calcu- 
lated. Indeed, AIT may be markedly distinct from the AII| 
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Fig. 3. The dipole (I = 1) asymptotic period spacing 
(AII|) in terms of the effective temperature for various 
stellar masses. Solid (dashed) lines correspond to stages 
before (after) the models reach their highest effective tem- 
perature (evolutionary knee). Also, the location of pulsat- 
ing PG 1159 stars with observed mean period spacings is 
shown. See Table 1 for details. 



predictions depending on the range of periods in which the 
average of the period spacing is performed (or observed). 
This is particularly true for the evolutionary stages cor- 
responding to the white dwarf regime, where, for a given 
stellar mass, AIT turns out to be about 1 s smaller than 
AITf when averages are taken on short period intervals. It 
is apparent that only in the case of long periods do the 
period spacings given by AII| resemble those predicted by 
AIT, i-e. the asymptotic conditions are nearly reached in 
this case. In view of this, we expect that for those pulsat- 
ing PG 1159 stars on the white dwarf cooling track, that 
usually exhibit short pulsation periods, the stellar mass in- 
ferred from AIT becomes substantially smaller than the 
stellar mass determined from AII^. We address this issue 
in the following section. On the other hand, for the stages 
before the evolutionary knee, the mean AIT values tend to 
be larger than the asymptotic AIT| ones. 



4. Mass determinations from the observed period 
spacings 

Here we employ the evolutionary models described pre- 
viously to infer the seismic mass of selected pulsating 
PG 1159 stars by comparing the asymptotic period spac- 
ing, AII^, and the average of the computed period spacings, 
AIT , with the observed mean period spacing, AIT 3 . These 
methods allow us to infer a value of the stellar mass as long 
as the effective temperature of the star is determined from 
spectroscopy or other method. Naturally one parameter, 
All, cannot determine two properties, Teff and logg. 

In Fig. [31 we show the evolution of AII| (for £ = 1) in 
terms of the effective temperature for the MA06 PG 1159 
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Fig. 4. The average of the computed period spacings for 
PG 1159 model sequences with different stellar masses in 
terms of the effective temperature. Each panel corresponds 
to a specific pulsating PG 1159. Also, the observed mean 
period spacings are shown. The top two panels correspond 
to evolutionary stages before the sequences reach the max- 
imum effective temperature, i.e., in the PNN stage instead 
of the DO stage. 



evolutionary models. The predictions corresponding to the 
evolutionary stages before the maximum effective temper- 
ature are indicated with solid lines, while the stages trac- 
ing the later evolution, hot white dwarf cooling branch, 
are denoted with dashed lines. In addition, the location 
of pulsating PG 1159 stars with the most recent deter- 
minations of the observed mean period spacings, All , 
are included in Fig. [3] — and also listed in the sixth col- 
umn of Table 1. Specifically, we include the observational 
data for PG 2131+066, PG 0122+200, PG 1707+427, 
RX J2117.1+3412, PG 1159-035, and NGC 1501. These 
pulsating stars are hot hydrogen-deficient, post-AGB stars 
for which the number of pulsation modes detected is 
high enough to infer an average of the period spac- 
ings. PG 2131+066, PG 0122+200, and PG 1707+427 are 
evolved PG 1159 stars on the hot white dwarf cooling 
branch and characterized by short pulsation periods (see 
last column in Table 1). By contrast, RX J2117.1+3412 
and NGC 1501, low-gravity and high -luminosity objects, 
pulsate with markedly longer periods. The pulsating 
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NGC 1501 belongs to the [WCE] class, Wolf Rayet-type 
central stars of planetary nebulae with emission lines and 
believed to be the progenitors of PG 1159 stars. 

From the asymptotic AII| diagram shown in Fig. [3l 
the stellar mass of the above mentioned pulsating PG 1159 
stars is assessed. The results are listed in the second column 
of Table 1. Note that the seismic masses as inferred from 
the use of the asymptotic approach differ by more than 10% 
from the spectroscopic masses (the spectroscopic masses 
are taken from MA06 and listed in the fifth column of 
Table 1 This difference is particularly true for the short- 
period variables such as PG 0122+200, PG 2131+066, and 
PG 1707+427, for which the seismic mass becomes about 
18% larger than the spectroscopic one. 

From the discussion in the previous section, we expect 
for our target stars smaller stellar masses when they are 
derived from the mean AIT. This is borne out by Fig 2J 
which displays AIT for t = 1 modes in terms of the ef- 
fective temperature for different stellar masses. Each panel 
corresponds to a specific star, discussed above. To derive 
the average of the period sp acing s in RX J2117. 1+3412 
and NGC 1501, we computed AIT for the high-luminosity 
(PNN) regime of the evolutionary sequence models, while 
for the remaining stars we compute values of AIT for the 
stages following the evolutionary knee for the PG 1159 
stars, i.e. the low-luminosity (DO) regime. Also for each 
star, the mean AIT is calculated by averaging the model 
period spacings over the corresponding period interval in 
which the periodicities are indeed observed. This is the rea- 
son for the fact that the curves are different in each panel. 
In the third column of Table 1 we list the resulting estima- 
tion of the stellar mass for the six stars. For those pulsat- 
ing PG 1159 characterized by short pulsation periods, the 
seismic masses as derived by this approach are appreciably 
lower — up to 0.06 Mq lower — than the values inferred by 
using the asymptoti c per iod spacing. As we mentioned, this 
is due to the mean AIT being typically 0.7—1.0 s smaller 
than the asymptotic AII^ when short periods are involved, 
i.e. for stages after the evolutionary knee. Thus, the discrep- 
ancy between seismic and spectroscopic masses is markedly 
alleviated when the average of the period spacings is used 
instead the asymptotic ones. Indeed, the seismic mass in 
this case becomes at most 6 % larger than the spectroscop- 
ically derived masses, except for the hot pulsating PG 1159 
star RX J2117. 1+3412, the spectroscopical mass of which 
is more than 20% higher than the asteroseismological mass. 

5. Discussion and conclusions 

This paper explores the systematic discrepancy between 
spectroscopical and asteroseismological masses of pulsat- 
ing PG 1159 stars. Our motivation is the result of Miller 
Bertolami & Althaus (2007) that such discrepancy should 
not be attributed to uncertainties in post-AGB tracks, but 
possibly to systematics in the asteroseismological mass de- 
termination methods. Recently, Quirion has pointed to one 
of us (M3B) that a possible opacity change resulting from 
the spread of He/C/O abundances in PG1159 stars could 
be a source of uncertainty in the location of the tracks. We 

1 We mention that in the case of PG 1159—035 and 
PG 0122+200 the uncertainty in the measured surface gravity 
translates into an uncertainty of ±0.1 Mq in the spectroscopic 
mass (MA06). 



addressed this issue by calculating sequences in which he- 
lium and carbon are changed in the whole envelope above 
the helium burning shell. We find that changing helium into 
carbon by an amount of 0.4 by mass shifts the track by only 
0.02 dex in effective temperature (being bluer if carbon is 
higher). This translates into a shift of only 0.005 and 0.015 
Mq for the spectroscopic mass near the 0.51 and 0.6 Mq 
tracks, respectively. Thus, the precise values of the He/C/O 
abundances do not seem to introduce appreciable changes 
in the masses derived by MA06. 

Specifically, we have concentrated on the seismic masses 
that result from a comparison of the observed period spac- 
ings with the usually adopted asymptotic period spacings 
(AII^) used in most mass determination of individual pul- 
sating PG 1159 and the better suited average of the com- 
puted period spacings (AIT). On the basis of full PG 1159 
evolutionary models that consider the evolutionary history 
of progenitor stars (MA06),and the ensuing internal chem- 
ical profile, we have shown that the derivation of the stellar 
mass using the asymptotic period spacing is not appropri- 
ate in the case of PG 1159 stars. In particular, we demon- 
strate that for those pulsating PG 1159 stars characterized 
by short pulsation periods, i. e., the pulsating PG 1159 
stars on the hot white dwarf regime (DOVs), the asymp- 
totic AIT| differs appreciably (by more than 1 s) from the 
mean AIT. Only in the case of variables with long periods 
(PNNVs), like the high - luminosity, log-gravity pulsating 
PG 1159 stars, do the g~ mode period spacings giv en by 
asymptotic AII^ resemble those predicted by mean AIT. 
This is expected because the asymptotic conditions are ap- 
proached in the limit of very high radial order k. 

For quantitative inferences, we have computed the seis- 
mic mass resulting from the employment of the asymptotic 
and the average of the computed period spacing for those 
pulsating PG 1159 which have a sufficiently large number 
of detected modes to infer an observed value of the mean 
period spacing. Our selected stars are listed in Table 1, to- 
gether with the stellar mass inferences. The employment 
of the asymptotic theory, in principle formally valid for 
chemically homogeneous stellar models at high radial in- 
dex k, overestimates the seismic mass by about 0.06 Mq in 
the case of very short period pulsating PG 1159 stars like 
PG 2131+066 and PG 0122+200. Because PG 1159 stars 
are expected to be chemically stratified, estimations of the 
stellar mass from mean AIT are more realistic than those 
inferred by means of asymptotic AITf. Indeed, stellar mass 
derived from the mean AIT are in good agreement with 
the mass values obtained from detailed period fittings. The 
discrepancy between asteroseismological and spectroscop- 
ical masses is markedly alleviated by the employment of 
the average of the computed period spacing instead of the 
asymptotic period spacings. 

In closing, a Fortran program to derive, from our evo- 
lutionary sequences, averages of the period spacing for 
arbitrary period intervals is available at our web site 
http: / /www.fcaglp.unlp.edu.ar/evolgroup. 
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Table 1. Stellar masses for selected pulsating PG 1159 stars as derived from the asymptotic and the average period 
spacings (second and third columns). The fourth column lists the stellar mass resulting from detailed period fittings, 
when available. The fifth columns displays the stellar mass as inferred from spectroscopy (from T e g and g values from 
Werner & Herwig 2006) . The sixth columns corresponds to the observed period spacings and the last column the observed 
range of periods for 1 = 1. All masses are in solar units. 
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References: (a) Corsico et al. (2007b); (b) Corsico et al. (2007a); (c) Corsico et al. (2007c); (d) Reed et al. (2000); (e) Fu et al. 
(2007); (f) Kawaler et al. (2004); (g) Vauclair et al. (2002); (h) Costa et al. (2007); (i) Kawaler et al. (1995); (j) Bond et al. 
(1996). (**) The two mass values result from considering that the star is either after or before the evolutionary knee. 
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